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The role of the multiorbital effects on the emergence of frustrated electronic orders on the tri¬ 
angular lattice at half hlling is investigated through an extended spinless fermion Hubbard model. 

By using two complementary approaches, unrestricted Hartree-Fock and exact diagonalizations, we 
unravel a very rich phase diagram controlled by the strength of both local and off-site Coulomb 
interactions and by the interorbital hopping anisotropy ratio t'/t. Three robust unconventional 
electronic phases, a pinball liquid, an inverse pinball liquid, and a large-unit-cell \/T2 x \/T2 droplet 
phase, are found to be generic in the triangular geometry, being controlled by the band structure 
parameters. The latter are also stabilized in the isotropic limit of our microscopic model, which 
recovers the standard SU(2) spinful extended single-band Hubbard model. 
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I. INTRODUCTION 

In analogy with frustrated spin systems, frustration of 
the charge interactions by the triangular-lattice geometry 
constitutes a favorable playground for the emergence of 
novel phases. In this context, the pinball liquid (PL) was 
proposed a decade agcP as an original Coulomb-induced 
charge-ordered metallic phase in the framework of the 
extended Hubbard model. This phase has no classical 
equivalent and takes advantage of quantum fluctuations 
in order to lift the massive degeneracy of interacting elec¬ 
trons on quarter-filled triangular lattices in the classical 
limit. While the existence of the PL is now a consoli¬ 
dated fact, as it has been demonstrated by several com¬ 
plementary techniques and approximation^^H^, such an 
interesting electronic phase has not been observed exper¬ 
imentally in the 6>-(BEDT-TTF)2X materials for which 
it was originally predicted. This can be attributed to 
the presence of other competing effects not considered 
in the idealized theoretical descriptions, most notably 
deviations from a perfectly isotropic triangular lattice, 
the interaction with the lattice degrees of freedonP and 
the presence of long-rang e tail s of the Coulomb repulsion 
beyond nearest neighbor which all favor insulating 
stripe-ordered states. 

There exist other classes of materials with layered tri¬ 
angular lattices and sizable electronic interactions which 
do present interesting charge-ordered phases whose origin 
is not fully understood. These include transition-metal 
oxides such as the layered cobaltates Naa,Co02, which 
exhibit complex electronic patterns which can be tuned 
by electron dopin^nHEl, and the triangular nickelates 
AgNiOj^^L^ and Ag2NiO^^, which show a threefold or¬ 
dered metallic phase with anomalous metallic proper¬ 
ties. Another interesting class is that of transition-metal 
dichalcogenides. In lT-TaS2, for example, the ordered 
state displays a marked Mott character induced b y charge 
modulations with a large periodicity of VTS x 


and various other periodicities are found in other com¬ 
pounds. What all these materials have in common is 
that electrons live in bands constructed from d atomic 
orbitals. Bridging the ideas of frustrated charge order 
from their initial domain of application (single-band, lay¬ 
ered organic conductors) to such d-electron compounds 
requires us to account for the presence of multiple bands 
and to move to electron densities not restricted to one- 
quarter filling. 

Multiband effects come in two different kinds. The 
first is related to interactions that are present within 
the d-electron manifold already at the atomic level, 
most notably the on-site Hund interactions acting on 
the magnetic degrees of freedonP^I. These are known 
to favor the emergence of high-spin states and have 
been shown to strongly enhance the effects of electronic 
correlation^^^"^^. Their ability to stabilize a PL phase 
with unconventional metallic properties at a filling of one 
electron per site has been explored very recently^. The 
second type of multiband effect, which can also lead to 
novel and original properties, is of kinetic origin and has 
to do with the microscopic form and symmetry of the 
interatomic electron transfers. One remarkable example 
is the hidden kagome symmetry and fiat bands which 
have been pointed out in the layered cobaltate^, and 
which could be related with the e xperim entally observed 
Kagome order in these compound d^^ l ^^ l 

The purpose of this paper is to study how such multior¬ 
bital kinetic effects influence the frustrated charge order 
on the triangular lattice. To this aim we employ an atom¬ 
istic tight-binding descriptiorP^ and set up a two-orbital 
extended Hubbard model where the band structure can 
be tuned systematically as a function of the microscopic 
transfer parameters. Since our main focus is the explo¬ 
ration of novel charge ordered phases, it is a good starting 
point to resort to a spinless electron description. This ap¬ 
proach has been shown to capture the correct ordering 
patterns realized in single-band modeld^^^ in the limit 
of strong local Coulomb repulsion, where the magnetic 
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energy scales are typically much smaller than the ones 
controlling charge ordering. It has also been successfully 
used to study charge ordering in the context of multi¬ 
band models for correlated oxideP^ in the ferromagnetic 
state. 

Our results, obtained here at half filling via unre¬ 
stricted Hartree-Fock (UHF) mean-field theory and exact 
diagonalization (ED) on small clusters, show the emer¬ 
gence of a rich panel of charge and orbitally ordered 
phases. Most notably, we find three robust unconven¬ 
tional charge-ordered phases, whose occurrence can be 
tuned by varying the multiorbital band structure param¬ 
eters. Two of these phases have peculiar properties since 
a fraction of the electrons forms a charge order with a 
threefold symmetry breaking, while the other fraction is 
free to move on the remaining sites of the lattice, form¬ 
ing a honeycomb structure. These phases are called the 
pinball liquid (originally found in quarter-filled lattices 
and obtained here at half filling) and the inverse pin¬ 
ball liquid (which can be viewed as the dual to the PL). 
The third unexpected phase found in this work is the 
large-unit-cell y/l2 x y/l2 droplet phase, also found in 
the isotropic limit where our model reduces to the spin¬ 
ful extended Hubbard model on the triangular lattice, 
where it was overlooked in previous studies. Such phases 
could be of relevance to a variety of triangular d-band 
electron systems such as the cobaltates, nickelates and 
dichalcogenides. 

This paper is organized as follows. The microscopic 
model and the two different methods of solution are de¬ 
scribed in Sec. II along with their respective advantages. 
The phase diagrams obtained from both methods upon 
varying the interaction parameters and the multiorbital 
band structure parameters are presented in Sec. Ill, to¬ 
gether with a detailed description and characterization 
of the different ordered phases. Our main results are 
summarized in Sec. IV. 


II. MODEL AND METHODS 

A. Spinless two-orbital extended Hubbard model 

To explore the influence of multiorbital effects on 
Coulomb-driven charge ordering on the triangular lat¬ 
tice we write the following spinless two-orbital extended 
Hubbard model: 

= - y] tjfdl^djr'+B.c. + U'^ni^rii^ 

{ij),TT' i 

+ y y] (nit + n4) (njt + nj4,), (1) 
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FIG. 1: Evolution of the density of states of the noninter¬ 
acting system with respect to the t'/t ratio. The sketches 
illustrate (top) the triangular lattice with two orbitals per 
site and sublattices A, B and C, corresponding to threefold 
charge ordering, and (bottom) the main symmetry points in 
the original and reduced Brillouin zone for three-sublattice 
ordering. 


the orientation of the bond (i,j) on the triangular lat¬ 
tice. The second term is an effective on-site Hubbard 
repulsion between electrons on different orbitals, with 
= dl^dir being the local electron density on orbital 
r. The third term is the Coulomb repulsion between elec¬ 
trons on neighboring sites, which constitutes the driving 
force for charge ordering. Note also that in the electro¬ 
static limit, if the presence (absence) of an electron at site 
i is interpreted as a spin up (down), we immediately see 
that this last term is responsible for the presence of geo¬ 
metrical frustration on the charge degrees of freedom on 
the triangular lattice, in analogy with magnetic systems. 

Besides its fundamental interest per 5e, such a two- 
band description applies to actual materials with com¬ 
plete t2g and partially filled eg shells. In AgNi02 {t^ge^g 
configuration, formal valence Ni^+), for example, the or¬ 
bitals are split between an eg doublet occupied by one 
electron, and a completely filled t2g triplet that can be 
neglected to a first approximation by virtue of the large 
crystal-field gap of ^ 2 e\^^. Since the twofold eg 
orbitals \d^z‘^_r2) and \dx2_y2) form a pseudospin-1/2, 
we label them, respectively, as | t) and | 1) and denote 
the corresponding Pauli matrices by r* with i = 1 , 2 , 3 . 
The transfer integrals along the lattice vectors ui^U2i and 
U3 = U2—U1 sketched in Fig. [T]can be exp ressed in terms 
of two independent parameters t and t' 


tui 


t 0 
0 f 


t2 h\ . _( t2 -h \ 
h u)^ "" V-^3 U ) 


which we study at a density of one electron per site. The 
first term describes d-electrons (creation and annihila¬ 
tion operators d\^ and dir) nioving with transfer inte¬ 
grals tJJ' which depend on both the orbital type and on 


witht2 = (t+ 3 t')/ 4 , ts = ^/ 3 {t—t')/A and ^4 = ( 3 t+t')/ 4 . 
In full generality we take f /t in the interval [— 1 , 1 ] (all 
values outside this interval amount to interchanging or¬ 
bitals a and b) and set t as the energy unit. Note that 
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for t' = t the kinetic term reduces to two independent in¬ 
stances of the triangular isotropic lattice, and the model 
becomes analogous to the single-band spinful extended 
Hubbard model. 

The effective on-site Hubbard repulsion term in Eq. 
0 describes the interaction between electrons on differ¬ 
ent orbitals. It obviously has direct relevance to the study 
of ferromagnetically ordered stateP^, in which case only 
one spin species is present and two electrons necessarily 
occupy two different orbitals. In a more general context, 
the present spinless model can be viewed as an approx¬ 
imation to tackle the strongly interacting (i.e. strong 
Hubbard repulsion and strong Hund coupling) limit of 
the two-band spinful mo del int roduced in Ref. [ 24 ]. Fol¬ 
lowing standard notation d^^ * ^^ *, the interaction energy of 
aligned-spin configurations on a site is U = U — 3 Jh^ 
where Jh is the Hund exchange coupling and U is the 
intraorbital Hubbard repulsion. Such configurations are 
favored (in other words. Hand’s coupling favors high- 
spin states) because all other configurations have ener¬ 
gies U — 2 Jh and higher. When Jh approaches t// 3 , 
the condition U — 2 Jh ^ U makes it is possible to re¬ 
strict the system to states with aligned spins since these 
become energetically more favorable, projecting out all 
other high-energy configurations. Spinless fermions con¬ 
stitute a reasonable approximation to this projection for 
those aspects of charge ordering which do not involve the 
magnetic degrees of freedom. Finally, a notable advan¬ 
tage of the spinless model is that it allows us to cross¬ 
check the mean-field results via the use of ED techniques 
which would be impossible if we considered the full or¬ 
bital and spin character in the model. 

We mention here that although related multiband 
models for charge and orbital ordering on the triangular 
lattice have been studied in recent years, how multior¬ 
bital kinetic terms affect frustrated charge ordering (in 
particular t he p inball liquid) remains an open question. 
Vernay et for example, studied the evolution of or¬ 
bital ordering as a function of t'/t in a spinful model via 
both mean-field and exact diagonalization, but they did 
not consider the charge orderiM induced by the intersite 
repulsion V. Uchigaito et al.^ performed a mean-field 
analysis of the effects of both the Hund and Jahn-Teller 
couplings on the electron ordering as a function of t'/t; 
however, because the repulsion V was not included in the 
model, no pinball-liquid phase was found at realistic val¬ 
ues of the Hubbard repulsion I/, which is at odds with 
the experimental observations in AgNi02. The question 
has also been addressed from an ab initio point of vievP^, 
including both local and nonlocal interaction effects, but 
without providing systematic studies as a function of the 
microscopic Hamiltonian parameters. Finally, multiband 
effects on charge ordering on the triangular lattice have 
been studied via both UHF and dynamical mean-field 
theory in Ref. [ 24 ] , but only the fully isotropic limit t' = t 
was explored. 


B. Methods 

To solve Eq. 0 we have used the UHF and ED meth¬ 
ods, which are both defined on the site basis. There are 
two types of clusters with periodic boundary conditions 
on the triangular geometry which respect all the symme¬ 
tries of the infinite lattice (translations and point-group 
symmetries). In terms of the basis vectors ui and U2^ 
they are defined by two vectors, 

Ti = lui + mu2 

T2 = -mui + (/ + m)u2, 

with / or m = 0 for regular clusters (considered here for 
the UHF computations) and / = m for tilted ones used 
for the ED. The regular clusters have a number of sites 
N = P, while it is for the tilted lattices. In 

order to allow all relevant symmetry breakings, we have 
considered lattices for which the number of sites is always 
a multiple of 12. 

In UHF, the interaction terms are decoupled as 

a Hartree term “ (^^t) ^ Fock 

term + dl^dji{d]^di^) - {dl^dj^){d]^di^). 

The sets (uir) and for i and j G [ 1 , A^], are com¬ 

puted from the wave functions in the single-electron ba¬ 
sis. A self-consistent loop is performed starting from dif¬ 
ferent initial trial states (homogeneous, random, charge 
and/or orbitally ordered) until a fixed point is reached. 
The result with the lowest energy is selected among the 
converged states. The method is free from local con¬ 
straints and ad hoe symmetrizations of the solution; 
that is, no particular form of the ground state is as¬ 
sumed. This allows us to obtain the most general or¬ 
dered states of the model in an unbiased way and has 
been proven v ery su ccessful in predicting novel phases in 
related model^^^^. 

Also due to the absence of constraints, the resulting 
mean-field Hamiltonian does not necessarily commute 
with and r^, which are, respectively, the total orbital 
pseudospin operator and its z component. We notice 
that the Fock terms of the Hubbard interaction can be 
recast in terms of the ladder pseudospin operators as 
allowing for the extraction of the aver¬ 
age orbital components (r^), (r^), and (r^). This is par¬ 
ticularly useful in order to characterize, in addition to the 
charge symmetry breaking, solutions having an orbital 
order such as the 120 phase already observed in the Hub¬ 
bard model at large uESMl. To avoid spurious solutions 
and to obtain smooth convergences, we have employed 
here a finite-temperature version of the UHF, checking 
for all our results that the low-temperature regime was 
reached and no thermal fluctuations remained. Typi¬ 
cally, this is always achieved for an inverse temperature 
(3 = l/ 50 t. Finally, we have performed a systematic size¬ 
scaling analysis for clusters up to 72 x 72 sites in order to 
identify transition lines in the thermodynamic limit. 

For the ED calculations, we have taken the largest 
available cluster fulfilling all the symmetry requirements 
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mentioned above, namely, the 12 -site tilted cluster (/ = 
m = 2, as depicted in Fig. Q. Since the Hamiltonian 
does not preserve the orbital flavor, one has to consider, 
for a fixed number of particles, all the orbital sectors. 
In addition, due to the structure of the hopping matrix, 
as soon as t' 7 ^ t, the three directions of the lattice be¬ 
come inequivalent hence breaking the Cqv point-group 
symmetry. The system can thus be block-diagonalized 
simultaneously for k points compatible with the remain¬ 
ing C 2 symmetries, which contain only two generators. 
Id and the rotation of an angle of tt around the 2 : 
axis. This is true for only points T and M; K and K/2 
are not compatible with such symmetries (see the corre¬ 
sponding Brillouin zone in Fig. [^. In a given k sector, 
the ground state (GS) can then be either symmetric (re¬ 
ferred to as A) or antisymmetric (Ap) against the C 2 
symmetries. This will be used in the interpretation of 
the ED phase diagram. 

For both UHF and ED we have identified the differ¬ 
ent phases via the following quantities: (i) the charge- 
charge correlations C{k) = (l/A/')('0o|p(—A:)p(A:)|' 0 o )5 
with p(k) being the Eourier transform of the total on¬ 
site density operator -h (ii) the orbital-orbital 
correlations T{k) = (I/A^)('0o|t^(—/c)t^(/c)|'0o), where 
is the ^ component of the pseudospin-orbital op¬ 
erator, and (hi) the average double occupation D = 
(I/A^)('0o| Einally, we also computed the 

spectral function 


system is a semimetal with two parabolic bands touch¬ 
ing at the F point and Dirac cones at the K points. 


B. Unrestricted Hartree-Fock 


a) 


> 


c) 






Ar{k,uj) =—lm{'tpo\dkT - 1 , ■ drl'^o) (2) 

TT UJ — rL ITj 

within the UHF approach to analyze the reconstruction 
of the band structures and Fermi surfaces. 


III. RESULTS 
A. Noninteracting limit 

Figure reports the evolution of the non-interacting 
density of states (DOS) as a function of t'/t. For t'/t = I, 
the model reduces to two independent instances of the tri¬ 
angular lattice. As soon as t' ^ t, however, the electronic 
dispersion separates into two nondegenerate bands. Cor¬ 
respondingly, the logarithmic singularity in the DOS of 
the triangular lattice splits into two peaks, and Dirac 
cones appear in the band structure at the K points. Al¬ 
though the DOS never vanishes at the Dirac points, due 
to the simultaneous presence of other bands at the same 
energy, these can be identified in Fig. by the kinks 
located between the two Van Hove singularities. A Lif- 
shitz transition occurs for t' 0.43t where the lowest of 
the two peaks in the DOS crosses the Fermi energy, and 
one of the two bands changes character from electron- 
like to holelike. The band structure eventually becomes 
particle-hole symmetric for t' = — t, in which case the 
Dirac point falls at the band center. In this case, the 


FIG. 2: (a) Phase diagram obtained from UHF at fixed 

t' — 0.5t. The phases denoted by HM and IPL are metallic, 
and all the other phases are insulating. IPL is charge ordered, 
PCOi and VT2 X \/T2 are charge and orbit ally ordered, and 
ODW and 120 are orbitally ordered, (b) Phase diagram at 
fixed U = I as a function of t'/t. HM, PL, IPL, and PCOm 
are metallic; all phases except HM are charge ordered, (c) 
Evolution of the three sublattice densities ua, nc and 
charge correlations C{M) along the constant U line shown in 
(a) [the scale of C{M) has been multiplied by 10]. For the 
\/T2 X \/T2 phase, tt-a, tib, nc are defined as the averages over 
the four nonequivalent sublattices in the new unit cell, (d) 
Sublattice densities as a function of t'/t for fixed U = I and 
U= I.I. 

The phase diagram obtained in the (U, U) plane for 
f/t = 0.5 from UHF is reported in Fig. [^a). We 
choose this representative value of f /t because the cor¬ 
responding noninteracting DOS qualitatively reproduces 
the main features calculated by DFT-LDA for AgNiO^p^. 
To explore all the possible regimes of the model, we also 
display in Fig. Hb) the phase diagram in the (t',U) 
plane at a fixed value of U = 1. We recall here that such 
a moderate U actually stands for a local Hubbard repul¬ 
sion U which is considerably larger than this value. As 
discussed after Eq. Q , U is an effective parameter which 
is related to the microscopic interaction parameters in 
the d-electron manifold via the equality U = U — SJh- 
Eor a choice of the Hund coupling Jh = 0.25U, for ex¬ 
ample, the value U = 1 corresponds to U = 4, a value 
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appropriate to moderately correlated materials. 

As we show below, a rich variety of phases is found 
and can be classified into charge ordered (at large V) 
and charge homogeneous (at low V). Among the latter, 



FIG. 3: Typical charge!- and orbital-density snapshots in 
the different charge-ordered phases realized in the model: (a) 
IPL, (b) PL, (c) X VT2, (d) PCOi, (e) PCOn, and (f) 
PCOiii. The radius of the disks is proportional to the charge 
density on the sites, and the light and dark fillings correspond 
to the partial orbital densities. 

several types have been identified in the model, whose 
charge and orbital density patterns are depicted in Fig. 

These are denoted as inverse pinball liquid (IPL), PL, 
pinball charge order (PCO), and the x ^/l 2 droplet 
phase. Supplementing the real-space snapshots. Fig. [^c) 
conveniently reports the evolution of the local charge den¬ 
sities as a function of V along the vertical line shown in 
Fig. I^a) corresponding to U = 1. An analogous scan 
is presented in Fig. id) for the t' dependence at fixed 
V = 1.1. 

a. Inverse pinball liquid. The IPL order develops 
in the whole region t'/t > 0.1 upon increasing the in¬ 
tersite repulsion V from the homogeneous metal phase 
(HM) for sufficiently low values of U [Figs, [^b) and 
[^c)]. This phase exhibits charge order with a three- 
sublattice structure (A, B, C). The charge density sepa¬ 
rates into charge-rich sites forming a honeycomb lattice 
{halls^ density = 1 + with 5 being the charge 

disproportionation) and charge-poor sites {pins^ density 
nc = 1 — 26) located on the remaining triangular sub¬ 
lattice (we conventionally label the sublattice densities 
in descending order, ua > tib > nc)- The charge dis¬ 
proportionation progressively increases with V towards 
the maximum allowed value 6 = 1/2 corresponding to 
a fully depleted charge-poor sublattice {nc = 0 and, 
correspondingly, ua = = 3/2). There is no orbital 

polarization in either the charge-rich or charge-poor sub¬ 
lattices. 

The IPL is metallic due to the presence of itinerant 
carriers (balls) on the charge-rich honeycomb network. 
Figure l^a) illustrates the Fermi surface (FS) obtained in 
this phase (V = 0.9, U = 1, t' = 0.5), which clearly shows 
the existence of holelike carriers around the K points of 


the original Brillouin zone (large hexagon), resulting from 
the folding of one of the noninteracting bands (shown as 
dashed lines). Small pockets can also be seen around the 
corners of the reduced Brillouin zone (small hexagons). 
These are remnants of the second band of the noninter¬ 
acting system, which at this value of t' crosses the re¬ 
duced Brillouin zone very close to its corners (denoted 
as K' points) and therefore folds into closed pockets of 
trigonal shape. The origin of the large hole FS and the 
small trigonal pockets can also be clearly seen by com¬ 
paring the spectral function A{k,uj) illustrated in Fig. 
0b) with the dispersion of the two noninteracting bands 
gashed lines). 

b. Pinball liquid. When t'/t < 0.1, the IPL |-|-0 
charge pattern shown in Fig. [^a) is replaced by the 
PL, with a 2-|-| pattern as illustrated in Fig. 0b). 
This phase is dual to the IPL in the sense that the roles 
of the charge-rich and charge-poor sublattices are inter¬ 
changed. We identify this phase with the original PL of 
Hotta and coworker^ because, from the point of view of 
the electronic densities {ua = 2, ub = nc = 1/2), it can 
be viewed as two realizations of the PL phase found at 
quarter filling, one per each orbital character (the PL at 
quarter filling has = 1, = nc = 1/4). Because of 

orbital-orbital interactions, however, the two realizations 
are not independent, and the present PL can occur for 
only small or moderate values of U ^5. This can be 
contrasted with the quarter-filled case, where a strong 
Hubbard term is required to stabilize the PlW^. 

The phase diagram in Fig. Ub) shows that the se¬ 
lection between the PL order and its dual IPL is en¬ 
tirely governed by kinetic effects. This can be under¬ 
stood by observing that the electrostatic energies of the 
PL and IPL are formally equal and do not depend on the 
sign of the charge disproportionation J, Ejpl = Epp = 
3H(1-(52). 

Similar to the IPL, metallic behavior arises in the PL 
from the motion of electrons living on the honeycomb net¬ 
work. At this value of t', the noninteracting FS (dashed 
line) is composed of a large electronlike surface around 
the F point and smaller holelike pockets around the K 
points. Contrary to the IPL shown previously, however, 
the overall shape of the FS is only weakly affected by 
charge ordering of the PL type because the original Fermi 
pockets do not cross the boundaries of the reduced Bril¬ 
louin zone [Fig. [^c) andj^d)]. 

e. Pinball eharge order. At large H, the system un¬ 
dergoes further charge ordering, stabilizing a 2-1-0 charge 
pattern termed pinball charge order [Figs. ^d)-|^f)]. 
Simple electrostatic arguments predict that the PCO is 
stabilized for H > L^/3, which nicely agrees with the nu¬ 
merical results at large U and V [see Fig. [^a)]. The 
threefold disproportionation splits the electronic disper¬ 
sion into three separate bands per orbital state. Count¬ 
ing the occupied states leaves us with a central half-filled 
band, which should lead, in principle, to a metallic be¬ 
havior. We find instead that the PCO phase found in the 
whole region t'/t > 0.1, denoted PCOi in Figs. [^a) and 
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FIG. 4: Fermi surfaces and spectral functions for the threefold metallic charge-ordered phases, obtained from Eq. § with 
a Lorentzian broadening 77 = 0.02 on 72x72 site lattices: (a) and (b) IPL {V = 0.9,U = l,t' = 0.5), (c) and (d) PL 
(V = 1.1,U — l^t' — 0), and (e) and (f) PCOm {V = 2,U — l^t' — —1). The solid lines in the top panels show the original 
and reduced Brillouin zone, whose symmetry points are labeled in Fig. ^ The dashed lines are the Fermi surfaces and band 
structures in the noninteracting limit. 
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FIG. 5: Spectral functions for the insulating charge-ordered phases: (a) PGOi (t'/t = 0.5, = 1 , E = 2.6), (b) PGOn 

{t'/t — 0, U = 1, V — 2), and (c) ^/V2 x ^/V2 droplet phase (F/t = 0.5, U = 1, V = 2.5). The dashed lines represent the band 
structures in the noninteracting limit. 


ib), is insulating [see the spectral function in Fig. Ha)]. 
This is ascribed to the presence of a spiral 120 orbital 
order on the singly occupied B sites, caused W the local 
interaction U, as seen in the snapshot in Fig. ^d). Other 
orbital orderings on the B sublattice are possible within 
the present 2-1-0 charge arrangement upon varying t'/t 
[Fig. [^b)]. These are PCOn for f /t ^0.1 and large V 
[uniformly polarized B sites, Fig. I^e)] and PCOm for 
^ fow ^ [unpolarized B sites, Fig. [^f)]. 

PC On is also an insulator because the ferro-orbital po¬ 
larization is sufficient to split the narrow half-filled band 
at the Fermi energy into two separate components [Fig. 
[^b)]. Note that this is the only phase in the whole phase 
diagram which supports a nonzero (t^).PCOiii instead 


has a metallic character. In the limit t' = —t, shown in 
Fig. B c), it is a semimetal with bands touching at the F 
and K points. 

If one interprets as usual the 120° ordering as the 
mean-field equivalent of a Mott insulating phase, it 
can be argued that the PCO is a neat example where 
charge ordering enhances the effect of Hubbard-type in¬ 
teractions, which happens because the reconstruction of 
the band structure leads to half-filled bands that are 
much narrower than in the homogeneous phase, as these 
rapidly shrink upon increasing V [Fig. |^a) and|^b)]. 
One notable example of this positive interplay between 
charge-ordering and the M ott mechanism is found in the 
dichalcogenide TaS 2 p^ (although the charge ordering 
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pattern there has a larger periodicity of 13 sites per cell, 
as discussed below). 

d. \/l2 X \/l2 droplet phase. For t'/t > 0.1, an ad¬ 
ditional phase emerges in the small-L^ regime, located 
between the two charge-ordered phases found at small 
and large V. This phase, which is stabilized with respect 
to the PCO phase by purely kinetic effects,!^ has a large 
periodicity with 12 sites in the unit cell. It is detected by 
the coexistence of peaks in the charge-charge correlation 
function at the K points (characterizing threefold order) 
plus all the M (whose combination describes a phase with 
fourfold symmetr}®) and K/2 points. The latter are 
absent in the phases with pure threefold symmetry. The 
M-point correlation is shown in Fig. [^c). The real-space 
snapshot in Fig. |^c) shows that the charges form discon¬ 
nected hexagonal droplets, composed of essentially dou¬ 
bly occupied sites. These are surrounded by sites which 
are either empty or almost empty. Correspondingly, this 
phase has a strongly insulating character due to the open¬ 
ing of a gap oc y at the Fermi energy [Fig. Igc)]. This 
phase also displays weak orbital order on the hexagons, 
signaled by nonzero orbital correlations at points M and 
KI2. 

We note that the emergence of a charge pattern with a 
high-order periodicity such as the one found here is not 
at all trivial given that the electron-electron repulsion is 
restricted to only nearest-neighbor sites. The region in 
parameter space where the \/l 2 x y/l2 droplet phase is 
stabilized actually includes the limit t' = t where Eq. 
0 becomes equivalent to the spinful extended Hubbard 
model, but it was overlooked in previous studie^^^^^®Il. 
Interestingly, the present 12-site period is very close to 
the 13-site star of David modulation found in the triangu¬ 
lar dichalcogenide compound IT-TaSj^^. The hexagonal 
droplets in Fig. [^c) are also similar to those recently 
predi cted theoretically in the kagome lattice at n = 1/3 
fillin^^^EH, but in that case they are a natural conse¬ 
quence of the larger unit cell of the underlying lattice 
(see also Refs. gSEl). 

Finally, we mention that the x ^/l 2 phase ob¬ 
tained in the narrow interval —1/3 < t'/t < —0.1 [see 
Fig. I^b)] has the doubly occupied hexagons replaced by 
empty sites and vice versa and can therefore be consid¬ 
ered the dual to the droplet phase described above. The 
two are separated by a phase with coexisting incommen¬ 
surate charge and orbital order, denoted C/ODW. 

e. Homogeneous orbitally ordered phases. Several 
phase transitions are also found within the charge ho¬ 
mogeneous region at low V. The system evolves upon 
increasing the effective local interaction from a paraor- 
bital metal at small U to a spiral 120 orbitally ordered 
insulating phase at large U (with the orbital moments 
arranged in planes perpendicular to the lattice as soon 
as t' 7 ^ t), going through an intermediate region with 
more complex orderings, denoted as ODW in Fig. |^a). 
Although a precise study of this intermediate region is 
beyond the scope of the present work, we would like to 
stress the following points: (i) Within our real-space un¬ 


restricted approach, we have recovered the two inter medi- 
ate mean-field phases of the spinful Hubbard modeP^^, 
namely an incommensurate orbital (spin) ordered phase 
whose wave vector evolves with U and a zigzag stripe 
compatible with the points k = (Ttt, 0). (ii) While previ¬ 
ous studies have looked for solutions breaking the symme¬ 
try at a single wave vector, we find solutions compatible 
with two or more coexisting k vectors, possibly a mean- 
field indication of a tendency to a structureless orbital 
liquid state, (iii) Away from the case of the isotropic 
Hubbard model, i.e., as soon as t' 7 ^ t , the zigzag stripe 
phase seems to disappear, and only the incommensurate 
regions with coexisting order survive, (iv) At t' = 0.5, 
the critical parameter Uc is decreased by 25% with re¬ 
spect to the isotropic case t' = t, in agreement with the 
corresponding reduction of the bandwidth. 

C. Exact diagonalization 

To ascertain if the variety of phases found at mean- 
field level in Fig. I^a) is robust against quantum fluctua¬ 
tions, we have performed a systematic study of the model 
equation 0 via Lanczos diagonalization on a 12-site clus¬ 
ter. Such a cluster is compatible with all the symmetries 
of the expected charge-ordered phases except for PCOi, 
which has a nine-site unit cell due to the presence of ad¬ 
ditional 120 orbital order. It is also compatible with the 
orbital orders found by UHF in the low-E region, except 
for the (Ttt, 0 ) zigzag stripes and the incommensurate 
stripes characteristic of intermediate U. 

Figure shows the phase diagram obtained in the 
(U^V) plane for f = 0.5t by combining the symmetry 
character of the ground state for the charge-charge [Fig. 
I^a)] and the orbital-orbital [Fig. [^b)] correlations at the 
corner of the Brillouin zone (K) and the average double 
occupation [Fig. [^c)]. Our ED results for small clusters 
confirm that the very rich physical picture created by the 
interplay between charge and orbital degrees of freedom 
persists even beyond the mean-field level. A large num¬ 
ber of domains with different symmetries are obtained, 
separated by first-order transitions, as displayed in each 
panel of Eig. these domains can be associated with all 
of the different phases found in Eig. [^a). 

a. Charge-ordered phases. The buildup of TC-point 
charge correlations in Eig. |^a), which is expected in all 
the charge-ordered phases described in Sec. HIB, shows 
remarkable agreement with the UHE charge-ordering 
transition lines reported in Eig. |^a) and asymptotically 
follows the analytical prediction W = U/3. We provide 
in Table 1^ some representative values of C{k) at points K 
and M for the different charge-ordered phases, as com¬ 
puted numerically by both UHE and ED and analyti¬ 
cally on ideal electrostatic patterns. As one can see, a 
quantitative agreement is found between the two numer¬ 
ical methods, which allows for direct identification of the 
mean-field phases in the ED results. Erom Eig. and 
Table [l| we associate the different charge-ordered phases 










FIG. 6: Phase diagrams obtained from ED on a 12-site cluster at fixed t' = 0.5t. (a) Charge-charge correlation C{K). (b) 
Orbital-orbital correlation T(K). (c) Expectation value of the double occupation D. The color maps are in arbitrary scale, 
increasing from dark to bright. On each panel are shown the ground-state symmetry sectors (translation and point group) and 
first-order transition lines (thick lines) reflecting symmetry breaking driven either by the charges, the orbitals or both. The 
symmetry sectors in the small dome at low V is respectively K/2 and as U increases. 




PCO 

IPL/PL 

VV2 X Vl2 

Analytical 

K 

1/3 

1/4 

1/4 


M 

0 

0 

1/9 

UHF 

K 

0.302 

0.121 

0.227 


M 

0 

0 

0.065 

ED 

K 

0.308 

0.157 

0.262 


M 

0.006 

0.055 

0.063 


TABLE I: Comparison of the charge-charge correlation 
C{k)/N for the high symmetry points K and M in the differ¬ 
ent charge-ordered phases, computed analytically on the ideal 
patterns, via UHE and ED. The values for PCO, IPL/PL, 
and X \/l2 correspond to t'/t = 0.5 and {U, V) of (6,4), 
(1,1.2), and (0,3), respectively. 


as follows: K/2 and K PCO, Ma^ PL, and Pa ^ 

v/l2 X v/l2. 

First, we identify both the iF/2 and K ground states 
with the broad PCO region of Fig. [^a). We attribute the 
additional transition seen here, which is absent in UHF, 
to the fact that the ninefold orbital order present in the 
PCOi phase is not compatible with the 12-site cluster, 
so that different orderings are stabilized instead. This 
may due to either the exact treatment of the correlations 
or the size of the lattice. Accordingly, this change of 
symmetry is not observable in C(iF), but it is clearly 
seen in T{K) in Fig. [^b). 

Second, by looking at the charge correlations in real 
space (not shown), it is possible to associate the Ma^ 
ground state with PL order. The PL is found here in¬ 
stead of the IPL expected from UHF because the con¬ 
sidered cluster is too small to capture the subtle kinetic 
effects which distinguish between these two phases. This 
is confirmed by the fact that the PL is also selected in 
the UHF solution when a 12-site cluster is considered, as 
we have checked. 


Third, the ^/l2 x \/T2 droplet phase found in the mean 
field can be associated with the F^ ground state at low U 
in Fig. This phase has a charge signature correspond¬ 
ing to a mixture of the high-symmetry k vectors 
and i^/2, with dominant weight on the first two, which 
agrees with the UHF and analytical results (see Table ^ . 
Also, it is in this phase that we find the strongest double 
occupancy [Fig. |^c)], which corresponds to the doubly 
occupied sites on the hexagons in Fig. ||:c). Our ED 
results confirm the finding that the y/l2 x vT2 droplet 
phase remains stable in the spinful extended Hubbard 
model (P = t). 

b. Orbitally ordered phases. The different phases found 
in UHF upon increasing U at low V also have their di¬ 
rect analogs in ED. We associate the spiral 120 phase 
obtained in UHF with the F^ phase at large U and small 
V\ here the double occupation is strongly suppressed, 
indicating a Mott insulating state with large thr eefold 
orbital correlations compatible with such an ordei^^Ell. 
We note that the critical value for the Mott transition, 
Uc ~ 9t, is smaller than the value Uc ~ 12t reported from 
analogous ED calculations in the isotropic casd^, which 
can be understood by observing that the bandwidth is re¬ 
duced by roughly 25% for the considered t' = 0.5t (Fi g§- 
We also observe a small dome in the middle of the Ma 
phase, which we tentatively associate with the ODW or¬ 
ders obtained by UHF on larger systems as suggested by 
the signatures in T{K). In this dome, two distinct sym¬ 
metry sectors corresponding to two different phases are 
encountered as U increases, K/2 and respectively. 


IV. CONCLUSIONS 

We have considered a minimal electronic model which 
describes the interplay between frustrated electron- 
electron interactions and multiorbital effects on the half- 
filled triangular lattice. Our results based on the com- 
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bination of a fully unrestricted Hartree-Fock method, 
which provides an accurate description of the multiband 
kinetic properties, and exact diagonalization on small 
clusters, which properly takes into account interorbital 
and inter site correlation effects, reveal a very rich phase 
diagram. A number of original charge-ordered and or- 
bitally ordered phases are displayed, whose occurrence 
can be tuned by varying the multiorbital band structure 
parameters and which could be relevant to a variety of d- 
electron systems on the triangular lattice. These include 
threefold metallic charge-ordered phases such as the pin¬ 
ball liquid, which was originally predicted to occur in 
quarter-filled lattices and is shown here to be stable at 
half filling, as well as its dual, the inverse pinball liquid. 
Both these phases could be relevant to AgNi02, where 
a robust threefol d char ge-ordered metal has been experi¬ 
mentally observecP^^^, which we therefore associate with 
a pinball state (see also Ref. [24]). An original insulating 
droplet phase with a large x \/l2 periodicity is also 
obtained here, which was overlooked in previous stud¬ 
ies of the extended Hubbard model and which could be 
closely related to the star of David charge-ordered phase 
of the triangular dichalcogenide compound IT-TaSj^. 
The insulating threefold PCO phase obtained here could 
also find a possible experimental realization, as suggested 
in triangular absorbate layerJ^. All these possible con¬ 
nections with the experiments emphasize the general na¬ 
ture of our study and should motivate further investiga¬ 
tions in stabilizing original charge orders driven by frus¬ 


trated electronic interactions and kinetic effects. 

Finally, we note that in this paper we have restricted 
our study to a perfectly stoichiometric case where the eg 
doublet is initially quarter filled, corresponding to an av¬ 
erage valence Ni^+ in AgNi02. It has been proposed, 
however, that the filling of the eg level s in nickelates 
might actually differ from this valu^^^^^, as extra elec¬ 
trons could be transferred from the oxygen 2 p orbitals 
provided that the electro static cost U = [/ — 3 J// on the 
Ni sites is sufficiently lovP^^^. The very observation of 
a pinball state in AgNi 02 indicates that the effective in¬ 
teraction U is low in this material [see Fig. 2(a)], and 
the existence of a “negative charge transfer” from the 
oxygen atoms is indeed compatible with the estimates 
of Ref. [16] , which indicate a reduced average valence of 
2.85+ per Ni. Studying the present model away from the 
perfectly quarter-filled configuration will certainly be of 
interest in view of these considerations. 
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